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ABSTRACT 

We compile a sample of 89 Seyfert galaxies with both [OIV] 25.89 /im line lumi- 
nosities observed by Spitzer IRS and X-ray spectra observed by XMM-Newton EPIC. 
Using [OIV] emission as proxy of AGN intrinsic luminosity, we find that although type 
2 AGNs have higher line equivalent width, the narrow Fe Ka line in Compton-Thin and 
Compton-Thick Seyfert 2 galaxies are 2.9^ 'g and 5.6^4 times weaker in terms of lu- 
minosity than Seyfert 1 galaxies respectively. This indicates different correction factors 
need to be applied for various types of AGNs before the narrow Fe Ka line luminosity 
could serve as intrinsic AGN luminosity indicator. We also find Seyfert 1 galaxies in 
our sample have on average marginally larger line width and higher line centroid en- 
ergy, suggesting contamination from highly ionized Fe line or broader line emission from 
much smaller radius, but this effect is too weak to explain the large difference in narrow 
Fe Ka line luminosity between type 1 and type 2 AGNs. This is the first observational 
evidence showing the narrow Fe Ka line emission in AGNs is anisotropic. The observed 
difference is consistent with theoretical calculations assuming a smoothly distributed 
obscuring torus, and could provide independent constraints on the dumpiness of the 
torus. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert - X-rays: 
galaxies - line: profiles 



1. Introduction 



Narrow Fe Ka emission line cores ar e common in the X- r ay spectra of both ty pe 1 and type 2 



Active Galactic Nuclei (AGNs) (e.g. see lBianchi et alJ 120091 ; lLaMassa et al.1 120091 for most recent 



studies). Such narrow lines have been traditionally associated with an origin in distant matter, 
especially the putative obscuring torus. Comprehensive and systematic studies of the narrow Fe 
Ka line in Seyfert 1 galaxies w ith t he Chandr a High Energy Grating (HEG) were presented by 
Yaqoob &; Padmanabhanl (j2004j ) and I Shu et al.1 (|2010l ). Shu et al. reported a mean line width of 
2060 km/s, supporting an origin in distant matter. Meanwhile, considering the observed larger 
scatter in line width from source to source, contributions from inner region (such as the outer part 
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of the accretion disk) could not be ruled out. Although the velocity widths in some of the sources 



were found consiste n t with those of the optical broad emission lines (also see lYaqoob et al.l 12001 



Bianchi et al.ll2003l ). iNandral (|2006l ) found no correlation between the Fe Ka core width and the 
BLR (specifically H/3) line width, eliminating the BLR as a general origin of the Fe Ka core (also 
see Shu et al. 2 010). The lack o f correlation between the equivalent width of narrow Fe Ka emission 
and EW(CIV) (jWu et all 12009 ) gives weight to this scenario. 



Type 2 AGNs are bel ieved viewed at inclination larger than type 1 AGN in the framework 
of unification model (e.g., lAntonuccil Il993l ). In such model, the line of sight of the observer to 
type 2 AGNs is obscured by an equatorial torus, which blocks the central accretion disk and the 
broad emission line region. The column density of the obscuration in type 2 AGNs spans from 
N H ~ 10 22 cm~ 2 to N H > 10 24 cm" 2 (Compton-Thick) (jBassani et al.lll999l ). The narrow Fe Ka 



emission lines often show higher equivalent w idth (EW) in type 2 AGNs, presumably due to the 
attenuation of the underlying continuum (e.g. lLaMassa et al.ll2009i ). Particularly, in type 2 AGNs 
with extreme X-ray obscuration (i.e. Compton-Thick), the narrow Fe Ka lines even appear as the 
most prominent featu re in the hard X-ray spectra and could have EW as large as several keV (e.g. 
Levenson et al .1 12006). As the narrow Fe Ka emission likely arises from the putative torus rather 
than the BLR, studies of the narrow core could place constraints on the unified model. 

Altho ugh the luminosity of the Fe Ka line has been suggested to be an indicator of intrinsic 
AGN flux (IPtak et al.ll2003l : iLevenson et al.ll2006l : lLaMassa et al.ll2009i ). it is still unclear whether 
the narrow Fe Ka line emission itself (i.e. the line flux but not line EW) depends on inclination. 
To answer this question, we present in this paper a first systematic comparison of the narrow 
Fe Ka lines in type 1 and type 2 AGNs. To perform such comparison, an intrinsic luminos- 
ity indicator of the central engine is required. Hard X-ray luminosity or reddenin g-corrected [O 



HI] were commonly used as isotropic luminosity indicat ors in various studies (e.g. iBassani et al 



1999 



Diamond-Stanic et al.ll2009i ; IMelendez et al.ll2008bl ). However, in the case of Compton-Thick 



sources which we are interested in particularly, hard X-ray luminosity is also strongly attenu- 
ated due to Compton scattering. The dust reddening to [O III] brings extra uncertainty to the 
measurement of intrinsic [ O III] luminosity. Furthermore, the [OIII] emission could be partially 
obscured in type 2 AGNs (jZhang et al.ll2008l ). making reddening correction procedure more dif- 
ficult. In this paper, we select to use [OIV] 25.89 fim line as an intrinsic luminosity indicator. 
Since [OIV] 25.89 fim has relatively higher ionization potential (54.9 eV), it's less affected by star 
formation. It is also si gnificantly less affected by ex ti nction than [OIII] (A v y 3 — 9 corresponds to 



i25.89Atm 



0.06-0.18: iGoulding k Alexanderll2009l ) . IMelendez et al.l (j2008bT ) have found that both 



[O IV] and [O III] luminosity correlate well with the very hard X - ray lu minosity (14-195 keV) as 
measured by the SWIFT/BAT. Meanwhile biamond-Stanic et al.l (|2009l ) reported that the [OIV] 
luminosity distributions are indistinguishable for obscured and unobscured Seyferts, while [OIII] 
luminosities are systematically smaller for obscured Seyferts. These confirm that [OIV] could serve 
as a reddening-free isotropic luminosity indicator. Also, [OIV] represents an improvement over the 
use of infrared continuum given the difficulty in isolating the AGN continuum from the host galaxy 



emission (jLutz et al.l 12004 ) 



2. Sample and Data Reduction 



2.1. Sample Selection 

[OIV] 25.89 fim emission lines have been detected in large number of AGNs by the In- 
frared Spectrograph (IRS) on board the Spitzer S pace Teles cope in the firs t Long- Low (LL1, 
19.5 ~ 38.0/im) IRS order. biamond-Stanic et all (hood ) and iMelendez etaD (ho08al ) presented 
two large AGN samples (respectively 89 and 103 Seyferts) with [OIV] measurements. We cross- 
correlated the two samples with 2XMM catalogue (the second comprehensive catalog ue of serendip 



itous X-ray sources from XMM-Newton) , which was released on 2007 August 22nd ([Watson et al 
2009). By excluding sources with less than 200 total EPIC counts in 2-12 keV, our final sample 



consists of 182 XMM observations of 89 Seyferts. We note that the final sample is never a ho- 
mogeneous one, since it relies on the availability of both archival Spitzer spectral data and XMM 
data of individual sources. However, most (if not all) of the Spitzer and XMM observations were 
planned independent of their OIV emission and Fe Ka line emission. Our study, which is focus on 
OIV and Fe Ka line emission, is thus free from the diverse selection effects. The cutoff in X-ray 
counts (> 200 counts, 2-12 keV) excludes weaker X-ray sources ( compared with OIV emission). 
However, this effect is also independent of Fe Ka line, since even for Compton-thick sources with Fe 
Ka EW of 1 keV, the emission line could only contribute 13% to 2-12 keV EPIC counts (assuming 
a pure reflection model pexrav with a photon index of 2 and a folded energy of 100 keV, plus a 
narrow gauss line at 6.4 keV whose EW is set to be 1 keV in xspec ). 



We gathered the redshifts and coordinates for each source from Simbad. X-ray classification of 
type 2 AGNs (to distinguish between Compton-thin and Compton-thick) requires not only proper 
spectral fitting to full band X-ray spectra, but also much more complex models than the ones 
adopted in this paper, including partial covering absorber, soft X-ray excess, scattering component, 
contamination from host galaxy, etc. Furthermore, low quality spectra could make the situation 
more difficult, and in many cases one need additional diagnostics, such as Fe Ka line EW and T 
ratio (F2-iokeV / F[oiii]> see Bassani et al. 1999). Thus a uniform summary of X-ray classifications 
is unavailable. Fortunately detailed X-ray studies for all but one Seyfert 2 galaxies in our sample 
have been published, and in this paper we choose to quote their X-ray classifications from the 
literature. 

The X-ray absorption column density of Seyfert 2 galaxy NGC 777 was not available from 
literature, and we classified it as Compton-Thin based on spectral fitting to XMM data. We di- 
vided the sample into 3 subsamples: Seyfert Is (including Sy 1.2-1. 5s, 33 sources), Compton-Thin 
Seyfert 2s (including Syl.8-1.9s, 35 sources) and Compton-Thick Seyfert 2s (with Nh > 10 cm~ 2 , 
21 sources). In panel A of Fig. Q]we plot the histogram distribution of [OIV] luminosity for three 
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subsamples. Excluding 6 upper limits of L[ iv]> we processed KS test (short for Kolmogorov- 
Smirnov test, which compares two samples under the null hypothesis that the samples are drawn 
from the same distribution Q) on [OIV] luminosity between subsamples. No significant differ- 
ence was found (the differences between Syls and Compton-Thin Sy2s, Syls and Compton-Thick 
Sy2s, Compton-Thin Sy2s and Compton-Thick Sy2s are at confidence levels of respectively 89.2%, 
14.6%, 85.5%, with corresponding D values of 0.29, 0.17, 0.32), ensuring we are comparing subsam- 
ples in same intrinsic luminosity range. Throughout this paper, we adopt cosmological parameters 
H = Whins' 1 Mpc~ l , fi m = 0.27, and fi A = 0.73. The sample is given in Table [TJ 



2.2. Data Reduction 

We used XMM-SAS (Science Analysis Software) version 9.0.0 and calibration files as of 2009 
August to reduce archival data. Each observation was processed using the pipeline "epchain" 
and "emchain". All the events were filtered to include only those with XMM-SAS quality flag 
FLAG==0 (#XMMEA_EM) for PN ( MOS1 and MOS2), and cleaned from flaring background 
with the task "espfilt". The task "espfilt" failed in 1 out of 9 exposures of Mrk 3, and this 
exposure (0009220301PNU002) was excluded from our further analysis. Another 4 observations of 
3C 273 in PrimeFull mode were eliminated due to serious pileup. Generally we extracted event files 
comprising single- and double-pixel events (PATTERN < 4) for PN, and single- to quadruple-pixel 
events (PATTERN < 12) for MOS1 and MOS2. However, in 90 exposures (Tabled]), moderate 
pileup was identified using the task "epatplot", t herefore onl y single-pixel events were extracted 



for these 90 exposures to reduce the pileup effect (Ballet 1999). 

The level of pileup was also estimated through PIMMS^I using observed spectra, and we found 
a maximum pileup fraction of 4% in PN detectors in our sample. Although such level of pileup 
could still alter observed spectral shape, the effect on the fitting to narrow Fe Ka line, which is the 
aim of this work, is negligible (also see §2.3). 

We defined the source region as a circle with radius of 40" centered at the source position, 
and the background regions as three circles with the same radii around the source region, being 
central symmetric as much as possible, and kept away from the CCD edges, the out-of-time events 
strips and other sources (as illustrated in Fig. [2]). When we were unable to avoid these disruptive 
features, the radii of the source or background regions were reduced down to 25" (see Table [I]). We 
generated the source and background spectra together with the appropriate redistribution matrix 
and ancillary response file from the source and background regions for each exposure, using the task 
"especget". Then for each source, to make the Fe Ka feature prominent as it may not appear in 



1 A statistically significant difference corresponds to a confidence level of 95% or higher 

2 htt p: / /heasarc. nasa. gov / docs /software /tools / pimms .html 
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individual observations, the PN spectra and the MOS1, MOS2 spectra from different observations 
were respectively combined, using the Ftools task "addspec", (see Section 3.2 for further discussion 
on the combination). The combined MOS1 and MOS2 spectra will be referred to as MOS hereafter. 

2.3. Spectral Fitting 

All spectra were grouped so that each channel contains at least one count, and spectral fitting 
were performed with C-statistic using Xspec version 12.4.0. As our major goal in this work is to 
measure the flux of the narrow Fe Ko line, we fitted the 5-10 keV spectra with a simple model 
uniformly. A Gaussian line was used to fit the narrow Fe Ka line, and the continuum were fitted 
with an absorbed powerlaw (zwabs*powerlaw+zgauss) . Some examples of our spectral fitting are 
given in Fig. [3j 

Whenever available, PN and MOS spectra were fitted simultaneously with same model param- 
eters, except for the powerlaw and Gaussian line normalizations, which were set free to account 
for the discrepancy among detector calibrations. Normalizations from PN spectra are presented 
in these cases. Pileup effects would be more significant in MOS detectors, which could alter the 
continuum spectral shape. For these MOS data suffering obvious pileup (see §2.2), continuum 
parameters (photon index and absorption) were also setting free. 

In Table [T] we list the best-fit parameter for the narrow Fe Ka lines. For 25 of the sources, 
the detection of narrow Fe Ka lines were insignificant (with F-statistic confidence level < 99%). In 
NGC7314 and I ZW 1, the centroid energies (rest-frame) of the detected gauss line are larger than 
6.4 keV at > 3a leve§, which seem to be due to ionized Fe lines, and extra gauss line at ~ 6.4 keV is 
statistically not required. For these 27 sources where the detection of the Fe Ka line was not signifi- 
cant, we provide upper limit to their narrow Fe Ka line flux by fitting with a gauss line with central 
energy fixed at 6.4 keV (rest-frame) and the a at 44 eV (the median a of the detected lines for the 
full sample, see panel B of Fig. [1]). Fixing the energy and a to the median values found for each 
source's corresponding sub-population does not alter our results in this paper. In several sources, 
extra gaussian lines were statistically required (with F-statistic confidence > 99%), most of these 
lines could be attributed to either broader Fe Ka line, ionized Fe Ka line or Fe K/3 line (see Table[2]). 

The spectral fitting to the narrow Fe Ka line could be significantly affected by continuum 
complexities. To quantify this effect, we generated 390 artificial spectra based on various contin- 
uum models plus narrow lines, and apply our uniform model (wabs*powerlaw+zgauss) to fit the 
artificial spectra between 5 and 10 keV. The models we adopted for simulations include A) partially 
absorbed powerlaw (pcfabs*powerlaw+zgauss model in Xspec, with covering factor varying from 
to 1, and Tin between 10 21 ~ 10 24 cm~ 2 ), and B) absorbed powerlaw plus reflection component 



3 3cr errors of centroid energies were calculated through xspec, but only 90% errors were presented in the tables. 
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(wabs*powerlaw+pexrav+zgauss in Xspec, with nn varying between 10 21 ~ 10 24 cm . and relative 
strength of the direct component to the reflected component varying from to 10). The line width 
and EW of the simulated narrow Fe Ko line were selected to cover the whole range of the final 
best-fit parameters in our sample (see Table [T]) . Through the simulations we find that by restricting 
the fitting to 5 - 10 keV range, we can recover the simulated narrow Fe Ko lines with systematical 
deviation < 2%. Adopting wider spectral range, such as 2 - 10 keV, would yield significant bias 
to the narrow Fe Kq line measurement (as large as 300%), due to the improper fitting to the 
underlying continuum. 



3. Statistical Distributions 

3.1. Centroid Energy, Line Width and Line Equivalent Width 

In panel C of Fig. [1] we first plot the centroid energy distribution for narrow Fe Ko lines in 
three subsamples. We find while the line centroid energies in type 2 AGNs are generally consistent 
with 6.4 keV, those in type 1 AGNs tend to be slightly higher. KS test indicates that the difference 
in the narrow Fe Ka line central energy distributions between type 1 AGNs and type 2 AGNs is 
statistically marginal (with a confidence level of 98.4% and D=0.385). We also derived a weighted 
mean centroid energy of 6.409±0.002 keV for type 1 AGNs, and 6.402±0.004, 6.401±0.002 keV for 
Compton-Thin and Compton-Thick Sy2s respectively. Higher centroid energies of the detected Fe 
Ka lines of Syls listed in Table 1 could be due to contamination from unresolved highly ionized 
Fe line or inner region (such as t he accretion disk) component which is Doppler broadened and 
affected due to general relativity ( Fabian et al. 1989H 4 ! and such components could be weaker in 



type 2 AGNs due to strong obscuration from the torus. 

We note that by studying large samples of type 1 AGNs with HETG observations, Yaqoob 
& Padmanabhan (2004) and Shu et al. (2010) have demonstrated that the centroid Fe Ka line 
energies in type 1 AGN strongly peak at 6.4 keV. This indicates that the contamination to the 
narrow Fe Ka line core is much weaker in HETG spectra due to its much better spectral resolution. 

We also compared the Fe Ka line widths between type 1 and type 2 sources. Due to the 
limited spectral resolution of XMM PN/MOS, many of the narrow lines are unresolved, and only 
upper limits were given (see Tabled]). The line width distributions are plotted in panel B of Fig. [H 
where can see slightly larger line width in type 1 AGNs compared with t ype 2 sources. Making use 



of ASURV (The Astronomy Survival Analysis; iFeigelson fc Nelsonlll985l ). which can be used in the 



presence of upper limit data, and will be used in this paper whenever upper limits exist in the data, 
we performed Peto-Prentice Generalized Wilcoxon test, and found that the line width distribution 
in type 1 AGNs differs from that of type 2 AGNs with a confidence level of 92%. This pattern is 



4 Modeling the broad component is beyond the scope of th i s pape r. Sy stematic studi e s of b road Fe Ka line emission 
from accretion disk are available in Brenneman fe Reynolds! (2009) and Nandra et al. ( 2007 ). 
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consistent with the centroid energy distribution, suggesting the narrow Fe Ka line in XMM spectra 
of type 1 AGNs is likely polluted by either highly ionized Fe line or broader emission from smaller 
radius. 



Consistent with previous studies (e.g. iBassani et al.lll999l ). we do find larger Fe Ka line equiv- 
alent width in type 2 AGNs (see panel D of Fig. [1]). Through Peto-Prentice Generalized Wilcoxon 
test ( ASURV) , we find the Fe Ka line equivalent width distribution in type 1 AGNs differs from that 
of Compton-thin and Compton-thick Sy2s with confidence level of 70% and >99.99% respectively. 



To examine the quoted X-ray classification, we reproduced the anti-correlation between the 
Fe Ka e quivalent width and the 2-10 keV luminosity normalized to the intrinsic AGN luminosity 
found by Bassani et al. (jl999l ) , replacing reddening corrected [OIII] luminosity with [OIV] luminos- 



ity as proxy of the intrinsic luminosity (see Fig. The 2-10 keV luminosities were calculated by 
fitting the 2-10 keV band spectra with two independent absorbed powerlaw and a gaussian line. We 
note such fitting could be unphysical to many sources, and is only valid to estimate the observed 
2-10 keV fluxes. In Fig. [H the concentration of Compton-Thick sources at the upper left cor- 
ner confirmed their classifications, including three Compton-Thick sources (NGC3982, NGC4501, 
NGC7674) with Fe Ka line non-detected (upper limits to Fe Ka line EW were plotted). 



Bassani et al.1 (119991) adopted a threshold of F2-iokeV /^\oni],corr < 1 to identify Compton- 
thick sources. lLaMassa et al.1 (|2010l ) reported an average ratio of Log(F yoiv]/^[Oiii],cnrr) = -0.82 for 



12/xm selected Sy2s. Thus log-Z^-iofcevV-^o/V] < 0-82 can also be used as a boundary for Compton- 
Thick sources in our sample. In Fig. [J] we plot log L 2 -iokeV I '^[oiv] = 1 an d log EWf b = 2.5 
(which include all Compton-thick sources) to delineate the Compton-thin/Compton-thick bound- 
aries. Note there are four Compton-Thin sources fall in the Compton-Thick region. They turn out 
to be "nearly" Compton-Thick with \ogn H > 23 (NGC1358: 23.4; NGC7479: 23.6; MCG-2-8-39: 
23.7; MRK273: 23.8. See references in Table. [T]). We conclude that the quoted X-ray classifications 
are generally consistent with the measured Fe Ka line EW and log £2— lOieeV /L\oiv] ■ Re-classifying 
Seyfert 2 galaxies based on Fig. [4] will not alter the results presented in this paper. 



3.2. Line Luminosity 

In the upper panel of Fig. [5] we plot the narrow Fe Ka line luminosity versus [OIV] line 
luminosity for our samples. We first perform Buckley- James linear regression on the log Lp e — 
logLyoiv] distribution for each subsample with ASURV. Six sources with only Ly OI y^ upper limit 
were shown in Fig. but excluded from statistical analysis since ASURV can not handle upper 
limits in both data sets. The best-fit slopes are: 

Syl-1.5: logL Fe = \ogL [OIV] x (0.89 ± 0.09) + 4.62 ± 0.40 

C-Thin Sy2: log L Fe = log L [OIV] x (0.92 ± 0.08) + 3.02 ± 0.36 
C-Thick Sy2: log L Fe = log L [OIV] x (0.73 ± 0.12) + 10.24 ± 0.45 
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We note that excluding the extremely faint Seyfert 1 M81 does not significantly alter the best-fit 
slope and other results presented in this paper. 

To explore the difference in Lyp^ between subsamples, we compare the residuals of logLi?e 
from the best-fit line of logLp e — logLtpiv] f° r Syls for each subsample. This is valid since 
the three subsamples show consistent Ly iv\ distributions (see §2.1 and Fig. |SJ). With Peto- 
Prentice Generalized Wilcoxon Test (ASURV), we found significant differences between Syls and 
Compton-thin Sy2s, and between Syls and Compton-thick Sy2s in the log-L^e residuals (with 
confidence level of 99.98% and > 99.99% respectively, see Fig. [6|). This implies the distribution 
of log Lp e — log Lipiv] a ^ so differs statistically significantly between Syls and Sy2s. The difference 
in \ogLp e residuals between Compton-Thin and Compton-Thick Sy2s is only marginal (with a 
significance level of 86%). By comparing the mean values of the log Lp e residuals in subsamples 
we find the narrow Fe Ka lines in Syls are 2.9^ 'g times stronger in terms of luminosity than 
that of Compton-Thin Sy2s, and 5.6^4 times stronger than that of Compton-Thick Sy2s. Simply 
comparing log Lp e j log LtQiv] between subsamples yields consistent results, that the narrow Fe Ka 
lines in Syls are 2.6^*6 ti mes stronger than that of Compton-Thin Sy2s, and 5.7^'g times stronger 
than that of Compton-Thick Sy2s. 

For comparison, we plot 6.4 keV monochromatic luminosity (Lp e /EWp e ) versus LtQiv] f° r 
three subsamples in the lower panel of Fig. [U Adopting the same methods used above, we find 
that the best fitting slopes for Syls, Compton-Thin Sy2s and Compton-Thick Sy2s are respectively 
1.03±0.07, 0.90±0.10 and 0.89±0.11, and the 6.4 keV monochromatic luminosity in Syls are 4.2+J| 
times stronger than that of Compton-Thin Sy2s, and 43.2^J < j ! 'y times that of Compton-Thick Sy2s. 
These decrements, if simply attributed to photo-electric absorption, require column densities of 
8.6±li x 10 23 cm 2 and 2.3^Q2 x 10 24 cm 2 respectively (calculated using a simple photo-electric 
absorbing model wabs*gauss in Xspec, with the line energy and a fixed at 6.4 keV and eV). 

To examine whether our fitting could be biased by spectral combination which could produce 
spurious results for variable sources, we repeat our spectral fit to individual observations. Discarding 
the observations in which narrow Fe Ka line is no longer detectable, we calculated the weighted 
mean line flux for each source. Fig. [7] compares the weighted mean line flux from individual 
observations with the line flux from the combined spectra, between which no obvious difference is 
found. 

4. Discussion 

By normalizing to [OIV], which serves as an isotropic indicator of AGN luminosity, we find 
that the narrow Fe Ka line in Compton-Thin and Compton-Thick Sy2s are 2.9^g'g and 5. 6^*4 
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times weaker than S yls respectively. We n ote that based on the X-ray reprocessing torus model of 



Krolik et al.l (|1994l ). iLevenson et al.l (|2006l ) performed numerical simulations to 7 Compton-Thick 
AGNs to estimate their intrinsic 2-10 keV luminosity from observed Fe Ka EW and luminosity. 
In their model, the obscuring torus is cylindrical symmetric, filled with constant density material, 
and has a square or rectangular cross-section allowing for unobstructed views over the half-opening 
angle. They obtained a typical value of Lp e / L intrinsic (2-iokev) = 2 x 1CP 3 for the 7 Compton- 
Thick AGNs. Such ratio is approximately 5 times smaller than the median value of our Syls with 
Fe Ka line detected (a median line EW of 96 eV could be converted to i^Fe/^-iofeeV = 9-7 x 10~ 3 
assuming a powerlaw with photon index T = 1.9). Taking our observed L/2 — lOkeV) for Syls 
presented in Table 1, we obtained a similar i^e/^b-iofceV = H-2 x 10~ 3 , 5.6 times larger than 
those for Compton-thick Sy2s presented by Levenson et al. These are remarkably consistent with 
our results that the observed Fe Ka emission in Compton-Thick AGNs is 5.6j^' 4 times weaker than 
type 1 AGNs. 

Meanwhile we also find marginally higher centroid energy and larger line width in type 1 
AGNs compared with type 2 AGNs. This indicates the narrow Fe Ka line in XMM spectra of type 
1 AGNs is likely polluted by highly ionized Fe line or broader component from smaller radius in 
type 1 AGNs, and such contamination is weaker in type 2 AGNs due to the torus obscuration. If 
highly ionized Fe line or broader component from smaller radius was responsible for the disparity in 
Fe Ka luminosity between Syls and Sy2s, the Syls with larger line widths and/or higher centroid 
energy should have systematically higher Lp e / 'Ly iv] ratios. In Fig. [6] we see that Syl galaxies 
with higher Fe Ka line centroid energy or larger line width ( compared with the median values of 
the Syls with Fe Ka line detected) do not have obviously larger Lp e /L\oiv\- KS tests show that 
the differences in logL^ e residuals in Fig. [6] between Syls with higher and lower centroid energy, 
and between Syls with larger and smaller line widths, are at confidence levels of 51% (D=0.31) 
and 91% (D=0.46) respectively, neither of which are statistically significant. Similar confidence 
levels of 51% (D=0.31) and 77% (D=0.38) were obtained if we simply compare Log(Lp e / Lujiv]) °f 
Syl subsamples. The average Log(Lp e / L^qiv]) 1S 0-18 ± 0.14 for Syls with higher centroid energy, 
—0.02 ± 0.09 for lower centroid energy, 0.16 ± 0.11 for Syls with larger line width, and 0.00 ± 0.12 
for smaller line width. This indicates that the contribution from either highly ionized Fe line or 
broader emission from much smaller radius is too small to explain the difference we found between 
type 1 and type 2 AGNs. 

We note that among the 26 Syls with narrow Fe Ka line detected in table 1, eight have both 
higher line centroid energy (E c ) and larger line width (<r), and eight show lower E c and smaller a 
( compared with the median values of Syl sample). Five Syls show higher E c but smaller <r, and 
the rest five have lower E c but larger a. Thus no significant correlation is found between E c and 
line width. Meanwhile, contamination from unresolved ionized Fe Ka emission could produce lines 
with higher E c and larger a, those lines with smaller E c and larger a are more likely due to Doppler 
broadening. 



The Fe Ka line emission could be determined by many factors, including the geometry, column 
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density, covering factor of the line-emitting gas, element abundances, and the orientation of the 
observer's line of sight relative to the line-emitting structure. In the standard unified model, the 
only difference between type 1 and type 2 AGNs is the orientation, thus is the only plausible factor 
to explain the different Fe Ka line emission. Observational selection effects could possibly favor 
smaller column density and covering factor of the torus in type 1 sources, but both effects predict 
weaker Fe Ka emission (contrarily, larger covering factor and column density of the torus could 
generally enhance the narrow Fe Ka line emission in type 1 sources). Meanwhile, we do not expect 
significantly higher Fe abundance in type 1 AGNs, not to mention that the Fe Ka line emission 
is ins ensitive to Fe abundance unless in optically-thin limit (N# < 10 23 cm -2 
2O10l ). 



see 



Yaqoob et al 



Yaqoob et al.l (J2010J) presented Monte-Ca rlo simulations on the pr oduction of Fe Ka emission 



in toroidal neutral X-ray reprocessor (also see lMurphy Yaqoobl 12003 ) . They clearly demonstrate 
that the strength of the line emission from a toroidal torus could dramatically depend on the viewing 



angle, as long as the column density of the reprocessor is significantly larger than 10 cm . For 
instance, for column density of 10 24 - 10 25 cm -2 , Fe Ka emission at face-on inclination could be 
around 3-30 times stronger than edge-on inclination (see Fig. 2 of lYaqoob et al.ll2010l ). This is 



mainly because of the larger optical depth to Fe Ka emission due to photo-electric absorption and 
Compton scattering at edge-on inclination. This scheme could easily produce the difference in the 
observed Fe Ka line emission between type 1 and type 2 AGNs. Clearly, this scheme requires the 
existence of Compton-thick torus in both type 1 and type 2 AGNs, although the line of sight is free 
of Compton-thick obscuration in type 1 and Compton-thin type 2 AGNs. 

However, if the line-emitting gas is clumpy instead of smoothly distributed, the inclination de- 
pend ency of the line emission could be reduced or even smeared out (also see the iNandra &: George 
19941 and iMiller et al.l [2009j) , dependent on the filling factor and covering factor of the clumpy 
blobs. Monte-Carlo simulations of the dust emission and radiative transfer in the clumpy torus 
have been performed to model infrared emission in various ty pes of AGNs, which could produce 



strong constraints on the dumpiness of the dusty torus (e.g. iNenkova et al.1 120081 ) . Our results 
could thus provide independent constraints to the distribution of the cold gas (which produces Fe 
Ka line emission, and its distribution may differ from that of dust) in a clumpy torus model once 
detailed Monte-Carlo simulation of Fe Ka line production in such model would be available. 

Our results indicate that the narrow Fe Ka emission in AGNs is anisotropic, with weaker 
emission (although higher EW) in obscured AGNs. Extra correction factors should be applied if 
one uses observed Fe Ka emission to estimate AGN's intrinsic luminosity. 



The work was supported by Chinese NSF through NSFC10773010/10825312, and the Knowl- 
edge Innovation Program of CAS (Grant No. KJCX2-YW-T05). 
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SI 
SI 
SI 



S2 
CT 
SI 



31 
26 



29 
3 

5 
4 



6 
7 
8 
9 
10 



30 
12 



Table 1 — Continued 



Name z ObsID Pileup;Radii log L2-l0fceV ^[OIV] ^ogLp e Ep e (keV) ap e (eV) EW (eV) Type Ref 



NGC 


4138 


0.002962 


0112551201 


1 


41.08 


38.92 


39.131°,;^ 


6.38l°;°i 


< 87.4 




S2 


8 


NGC 


4151 


0.003319 


0112310101 
0112830201 
0112830501 
0143500101 
0143500201 
0143500301 
0112190201 


1 

1111,1X12; 

ml,m2; 
ml,m2; 
•> 


42.33 


40.70 


40.6lt8;§§ 


6-401°;°° 


48.5+ 2 ;t 


11711 


SI 




NGC 


4168 


0.007388 


0112550501 


■> 


40.12 


39.23 


<38.98 


6.40* 


44* 


<901 


S2 


8 


NGC 


4235 


0.008039 


0204650201 


7 


41.65 


39.79 


Ar . _, ,_i_n 1 ^ 

40.14t°i| 


6 41 i u - u ^ 
-0.03 


o± - z -51.0 


. _|_1 1 ^ 

284 l98 


SI 




NGC 


4258 


0.001494 


0110920101 
0203270201 
0059140101 
0059140201 
0059140401 
0059140901 


) 

;PN:25 

t 
f 

: 
i 


40.54 


38.57 


38.241°;?° 


6 40+ ' 03 


< 98.8 


34+20 


S2 


8 


NGC 


4378 


0.008536 


0301650801 


i 


<40.54 


39.47 


<39.45 


6.40* 


44* 


<1110 


S2 


8 


NGC 


4388 


0.008419 


0110930701 


i 


42.55 


41.61 


. ,-.,-,-1-0 11 

41.00!;;^ 


_ . o 4-0 02 

6 - 43 Iao2 


49 fi+27.6 


156l« 


S2 


13 


NGC 


4395 


0.001064 


0112521901 
0112522701 
0142830101 
0200340101 


i 
i 
i 

; 


40.18 


38.02 


38.201°;^ 


f\ 37 ' u - uz 
D - J '-0.03 


— . ~^3X 7 
74 Q ' 


,-,-,-1-94 
S7--23 


S2 


8 


NGC 


4472 


0.003326 


0112550601 
0200130101 


;PN:25 

7 


39.91 


<39.21 


<38.17 


6.40* 


44* 


<229 


S2 


8 


NGC 


4477 


0.004520 


0112552101 


t 


39.64 


38.88 


<38.81 


6.40* 


44* 


<1436 


S2 


8 


NGC 


4501 


0.007609 


0106060601 
0112550801 


1 

;M1M2PN:25 


40.13 


39.71 


<38.89 


6.40* 


44* 


<871 


CT 


14 


NGC 


4507 


0.011801 


0006220201 




42.61 


41.01 


41 1Q+0.05 


6 3Q+ - 01 


48 2 +12,2 




S2 


15 


NGC 


4565 


0.004103 


0112550301 


;M1M2PN:25 


39.96 


38.89 


<38.65 


6.40* 


44* 


<386 


S2 


8 


NGC 


4579 


0.005067 


0112840101 


1 


41.39 


39.21 


on 57 +0.17 


fi oq+0.04 


n 7 +86.3 


137+^ 


S2 


16 


NGC 


4593 


0.009000 


0109970101 


m2; 


42.65 


40.38 


40.9lt°;i; 


6 42+ ' 01 
O-^_0.02 


o i+ 41 - 6 

u - 1 -0.2 




SI 




NGC 


4594 


0.003416 


0084030101 


1 


40.63 


38.83 


<38.99 


6.40* 


44* 


<243 


S2 


17 


NGC 


4639 


0.003395 


0112551001 


1 


40.18 


38.59 


<39.08 


6.40* 


44* 


<871 


SI 




NGC 


4698 


0.003342 


0112551101 




39.19 


38.70 


<38.09 


6.40* 


44* 


<834 


S2 


8 


NGC 


4939 


0.010374 


0032141201 


1 


42.02 


41.01 


<40.63 


6.40* 


44* 


<248 


S2 


18 



Table 1 — Continued 



Name 



ObsID 



Pileup; Radii 



log L Fe 



NGC 


4945 


r\ rim oto 

0.U01878 


(JllzolUoUl 


;JVLlM2JrJN :zo 


40.19 


39.37 








0204870101 


;JVLlM2JrJN :zo 






NGC 


4968 


0.009863 


0200660201 




40.86 


40.81 








0002940101 


; 






NGC 


5005 


0.003156 


0110930501 




39.98 


38.64 


NGC 


5033 


0.002919 


0094360501 




40.97 


39.48 


NGC 


5128 


0.001825 


0093650201 


ml,m2,pn; 


41.94 


39.86 








0093650301 


m2,pn; 






NGC 


5194 


0.001544 


0112840201 


i 


39.39 


39.11 


NGC 


5256 


0.027863 


0055990501 


i 


41.81 


42.04 


NGC 


5273 


0.003549 


0112551701 




41.35 


39.01 


NGC 


5506 


0.006181 


0013140101 




42.80 


41.27 








0013140201 


ml,m2; 












0201830201 


1 












0201830301 














0201830401 














0201830501 








NGC 


5548 


0.017175 


0109960101 


m2; 


43.45 


41.00 








0089960301 


i 












0089960401 


i 






NGC 


5643 


0.003999 


0140950101 


;M1M2PN:25 


40.48 


40.46 


NGC 


6240 


0.024480 


0101640101 


i 


42.51 


41.82 








0101640601 


i 












0147420201 


i 












0147420401 


i 












0147420501 


i 












0147420601 


i 






NGC 


7172 


0.008683 


0147920601 




42.58 


40.91 








0202860101 


i 






NGC 


7213 


0.005839 


0111810101 


m2; 


42.20 


39.20 


NGC 


7314 


0.004763 


0111790101 


ml,m2; 


42.23 


40.39 


NGC 


7469 


0.016317 


0112170101 


ml,m2; 


43.13 


41.34 








0112170301 


ml,m2; 






NGC 


7479 


0.007942 


0025541001 




<40.55 


40.57 


NGC 


7582 


0.005254 


0112310201 


1 


41.27 


41.13 



E Fe (keV) a Fe {eV) EW (cV) Type 



Ref 



,+0.04 
2 -0.03 



40.121°;" 



<38.63 
39-3llg 
39.85_ ; 07 



40 57+ ' 17 
,u,o '-0 .22 

39.681°-^ 
40 56 +o: ° 8 



41.32 



+0.06 
0.05 



39.70l°?o 
41 26+ '" 7 



40.13 



J -0.09 

+0.07 
0.07 

<40.17 
41.24l»;°5 

39.8ll°i° 

40.0518:°* 



6.40 

6.42 

6.40 
6.42 
6.42 

6.43 
6.46 
6.43 
6.40 



+0.00 



+0.02 

0.03 



+0.02 
0.03 
+0.01 
0.01 



+0.02 
0.02 
+0.03 
0.03 
+0.03 
0.03 
+0.01 
0.01 



6.41 



+0.01 



6.41 
6.41 



+0.01 
0.02 

+0.02 
0.01 



,+0.01 
-0.01 



6.42 
6.40* 



+0.01 
0.01 



6.41 



+0.02 

0.02 



6.45 
6.41 



+0.04 
0.04 
+0.01 
0.01 



< 4.4 

< 65.2 

44* 

< 67.4 

30 7+ 17 

' -30.7 



< 48.9 

< 40.6 



65.4 



+ 19.2 

20.4 



8.3_ 
44* 
64.1 



+36.8 
8.3 



+21.9 

29.1 



78.6lf 8 ;2 

1-13.4 
-29.1 



696lf 2 
1252l 3 gg 



<431 

206lg 

50l? 



7031-1 
418l?8 6 5 

9-1 o+83 

40l 7 



1139 
346 



+275 
240 
+57 
53 



i+13 
5 -16 

-13 
-10 

<64 

3 + 18 
-18 



74^ 



99^ 



29. r 



1016 

343 



+604 



CT 

CT 

S2 
SI 
S2 

CT 
CT 
SI 
S2 



SI 



CT 
CT 



+31 

33 



S2 

SI 
S2 
SI 

S2 
CT 



25 
19 
20 

21 

25 

22 



6 
32 



18 



23 



Table 1 — Continued 



Name 



ObsID 



Pilcup;R,adii logL2- 



logL 



[OIV] 



;L Fe E Fe {keV) a Fe (eV) EW (eV) Type Rcf 



0204610101 

NGC 7590 0.005255 0112310201 
0204610101 

NGC 7603 0.029524 0066950301 
0066950401 

NGC 7674 0.028924 0200660101 
UGC 12138 0.024974 0103860301 ;M1:25,M2:35 



40.12 39.62 
43.59 41.15 



42.20 
43.13 



41.95 
41.20 



<38.98 

<42.06 

<41.20 
<41.78 



;.40* 
1.40* 

;.40* 
;.40* 



44* 

44* 

44* 
44* 



<632 

<302 

<531 
<515 



CT 
SI 

CT 

S2 



21 



33 
19 



Note. — Column(l): Name of object. Column(2): redshift. Column(3): XMM-Newton observation ID. Column(4): before the semicolon: which of the 3 
EPIC detectors' exposures are suffering significant pileup. After the semicolon: in which of the 3 EPIC detectors' exposures, source region radius are reduced, 
instead of being the default 40 arcsec. Ml and M2 are abbreviations for MOS1 and MOS2. The number after the colon represents the extraction region radius 
in arcsec. Column(5-6): 2-10 keV and [OIV] luminosities in erg/s. Column(7-10): luminosity, central energy, a and equivalent width of the narrow Fe Ka line 
in rest-frame with 90% errors. The upper limits to a are at 90% confidence level. 3ct upper limits of luminosity and equivalent width are given when the line 
is not detected. * Denotes fixed quantities. Column(ll): Seyfert type. SI: Seyfert 1; S2: Compton-Thin Seyfert 2; CT: Compton-Thick Sy2. Column(12): 
References for the X-ray classification of Sy2 galaxies (Compton-thin or Compton-thick). 



Re fere nces 
Il998l: f7lAwaki et alj |2008|: 



(12 
(18 
|20q3 



this wor k: (2lTreister et all |2009|; OlChurazov et al.l |2003|: (4lRisaliti et al.l I2009I: (5lCardamone et alj I2007I: (6lMaiolino et"al 



_ .20081: (8lAkvL._ 
Shu et al.l 2007; (13TBeckmann et alj|2004 (14' Brightman fe Nandral 



Ge^rean topoulosl 120091 : (9lBre^mema*n fe Re^m^ld*sl~ l20*0*£ ; (lo]^omo^^^^_ 

"2008I: (15 lMatt et aTll2004l : (16toewangan et al 



Noguchi et al.ll2009l: ( 19 Gallo et al.ll2006l: (20lRothschild et al. 20061: (2llFukazawa etal 



24lBassani et al.l 19991: f25lGuainazzi et alj |2005|: (26lShu et al 



2008; 



2001 

(27lBraito et alj 2i M I 



(30 lKomossa fe Finklll997l ; (3l 1lwasawa et al.ll200ll ; f32 lVignati et al.lll999l; (33 lMalaguti et al.lll998l: 



22lShino 
iol (28 j 



inozaki et al. 



2001 



2006 



Balestra et al.l 



2004 



(23 



2005 



(11 
ill 



Madejsk^^alJ 



Pellegrini et al. 



2006 



2003 



Branduardi-Raymont et al. 
(29 lGreenhill et alj |200S| ; 
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Table 2. Additional Gauss lines 



Name Energy a Flux References 

(fceV) (eV) (lO-^ftotons/cm- 2 /;?- 1 ) 



3C 120 


6 96+ 03 

D - MD -0.03 


< 72 


9l 5 3 


Ballantvne et al. (20041 


Circinus 


7 04+ - 01 
'• u ^-o.oi 


39± 7 2 


41±* 


Molendi et al. (2003) 




7 47 +0.01 
'•*'-0.01 


< 31 


lit? 






6.69±°:gi 


67+ 14 

°'-14 


26+3 




I ZW 1 


6.83±g;15 


01 r) + 128 

Jlz -115 


8+? 


Gallo et al. (2007) 


MCG-6-30-15 


5 O7+0.09 
°- a '-0.05 


o 10 +98 
Jlu -102 


241- 


Fabian & Vauehan f 20031 




fi 8Q+ - 01 
D.89_ Q 02 


< 27 


e±i 


MRK 335 


7 n9 +0.03 
'■ uz -0.04 


< 58 


°-2 


Gondoin et al. (20021 


NGC 1068 


6 69+ 01 


63± 14 


33+3 


Pounds & Vaughan f 20061 




7 oi +0 - 02 

'■ ul -0.02 


68l a 1 


15±| 


NGC 1275 


6 69+ 01 
°- D9 -o.oi 

fi Q8+ 03 
b - 9 8_ .03 


< 16 




Churazov et al. f 20031 




< 58 






NGC 1365 


r 70+0. 09 


435±|| 


01 -24 


Risaliti et al. f 20091 


NGC 3516 


6 42+ 02 
°'^ z -0.01 


1271^ 


36+ 8 


Turner et al. (20021 


NGC 3783 


7 n o+0.02 
' ■ uo -0.02 


67 +22 
D '-25 


16l 3 


Reeves et al. (20041 


NGC 4151 


™8t$ 


< 65 




Schurch et al. (20031 


NGC 4945 


6.68t°;« 3 3 


218±?| 


18+ 2 


Schurch et al. (20021 


NGC 5506 


"■ OJ -0.06 


273l 3 J 




Matt et al. (20011 




6 99+ 03 
"■ a; '-0.02 


< 73 


11+ 6 




NGC 6240 


6.67±S:?,i 


< 54 


y+2 


Netzer et al. (20051 


NGC 7314 


6.68i:S 


3601& 


4,5+26 


Yaqoob et al. (20031 



Note. — Column(l): Name. Column(2-4): Centroid energy, <r and flux of the line. The errors 
and upper limits are in 90% confidence level. Column(5): References in which the detections of the 
additional lines were also reported. 
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6.55 
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10 
8 
6 
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2 
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Fig. 1. — The histograms of [OIV] line luminosities and centroid narrow Fe Ka line energies, a and 
EW. In panel B, C and D, we plot the median values and errorbars of the sources with narrow Fe 
Ka line detected. 
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Fig. 2. — Example of source and background regions used to extract spectra (NGC 2992, ObsID 
0147920301, PN detector), with single-pixel events only. The colorbar refers to the photon counts 
in each pixel. 
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£ 4x10" 
8 2x1 Q a 



Energy (keV) 



Energy (keV) 



A 




Circinus C-Thick Sy2 








A 


5 




Energy (keV) 






Energy (keV) 
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Fig. 3. — Sample spectra of our sample, including two Seyfert 1-1.5 galaxies, two Compton-thin 
Seyfert 2 galaxies, and two Compton-thick Seyfert 2 galaxies. The black and red data points refer 
to PN and MOS spectra respectively. All the model components but the narrow Fe Ka line are 
over plotted to show the Fe Ka line in the data to model ratio plot. The upper three sources have 
Fe Ka line detected, and the lower three do not. 
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Seyfert 1-1.5 
Compton-Thin Seyfert 2 
Compton-Thick Seyfert 2 



-1.0 -0.5 0.0 



0.5 1.0 1.5 

logL 2 _ l0keV / Lpj^ 



2.0 



2.5 



3.0 3.5 



Fig. 4. — Narrow Fe Ka line equivalent width versus / 'Lyoivy The yellow lines can be 
used as plausible boundaries for Compton-Thick sources. 
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T 




log LpiV] (erg/ s) 

Fig. 5. — Upper panel: Narrow Fe Ko line luminosity versus L^qiv]- For the sources with narrow 
Fe Ka lines detected, 90% errors of the luminosities are plotted. When the narrow Fe Ka line is not 
detected, 3cr upper limit of the luminosity is given. Lower panel: 6.4 keV monochromatic luminosity 
versus £[ozv]- In both panels, the best-fit lines for Syls, Compton-Thin Sy2s and Compton-Thick 
Sy2s are shown as blue, green and red lines. 



-27- 



o 

e 



o 



1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 
-2.0 
1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 

■ 2 -1 



O Syls, a > Median(a) • other Syls 




■ Compton—Thin Sy2 
a Compton— Thick Sy2 



? 




O Syls, E>Median(E) • other Syls 



37 38 39 40 41 42 

l °SL [0 ^(erg/s) 



43 



Fig. 6. — The deviation of log Lp e from the best-fit log Li? e -log L [iv] nne of Seyfert 1 galaxies, 
where we clearly see weaker narrow Fe Ka emission line in Seyfert 2 galaxies. Best-fit \ogLp e — 
log Lq[iv] correlations for Compton-thin and Compton-thick Sy2s are over-plotted. We also divide 
Syls by the Fe Ka line width and centroid energy using the median value of the Syls with Fe Ka 
line detected as a boundary, which demonstrate that Syls with larger Fe Ka line width or centroid 
energy do not show obviously stronger line emission. 
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Fig. 7. — For sources with multiple XMM observations, spectral fitting to composite spectra yield 
narrow Fe Ka line flux consistent with the weighted mean of individual exposures. The line repre- 
sents a slope of unity. 



